Apparent low-energy scale invariance in two-dimensional Fermi gases 
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Recent experiments on a 2d Fermi gas find an undamped breathing mode at twice the trap 
frequency over a wide range of parameters. To understand this seemingly scale-invariant behavior 
in a system with a scale, we derive two exact results valid across the entire BCS-BEC crossover 
at all temperatures. First, we relate the shift of the mode frequency from its scale-invariant value 
to 7^ = (1 + 2/d)P — p(dP/dp) s in d dimensions. Next, we relate 7^ to dissipation via a new 
low-energy bulk viscosity sum rule. We argue that 2d is special, with its logarithmic dependence 
of the interaction on density, and thus 72 is small in both the BCS and BEC regimes, even though 
P — 2e/d, sensitive to the dimer binding energy that breaks scale invariance, is not. 
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Systems exhibiting scaling symmetry or conformal in- 
variance are very special. In all laboratory realizations, 
one needs to tune one or more physical parameters (tem- 
perature, chemical potential, coupling) to observe scale- 
invariant behavior, for instance, in the vicinity of a quan- 
tum critical point [1]. Another example is provided 
by strongly interacting Fermi gases in three spatial di- 
mensions (3d), which display remarkable scale invari- 
ance properties at unitarity, where the s-wave scatter- 
ing length diverges by tuning to the Feshbach resonance. 
This is manifested in universal thermodynamics [2], the 
vanishing of the d.c. bulk viscosity [3], and the entire 
bulk viscosity spectral function £(w,T) [4] at unitarity. 
There may also be tantalizing connections between the 
ratio of the shear viscosity 77 to the entropy density s 
of the unitary Fermi gas [5, 6] and the bound for rj/s 
conjectured on the basis of gauge/gravity duality [7]. 

For the unitary gas in a 3d isotropic harmonic trap, 
scale invariance manifests itself most dramatically as an 
undamped monopole breathing mode oscillating at twice 
the trap frequency u>o independent of temperature [8, 9]. 
This mode corresponds to an isotropic dilation of the gas 
wherein the coordinates in the many-body wavefunction 
are scaled cx cos(wt). Scale invariance implies that this 
wavefunction is an exact eigenstate of the Hamiltonian 
and oscillates at a frequency 2wo without damping [8]. 

In a recent experiment [10], collective modes in a two- 
dimensional (2c?) Fermi gas were measured over a broad 
range of temperatures and interaction strengths. Re- 
markably, the breathing mode was found to oscillate 
without any observable damping at ~ 2ujq for 0.37 < 
T/T F < 0.9 and < ln(fc F a 2 ) < 500, where a 2 is the 
2d scattering length. This observation is extremely sur- 
prising, given that there is no a priori reason to expect 
scale-invariant behavior in a system which has a scale, 
namely, the dimer binding energy in 2c?. 

Our goal is to understand why the 2d Fermi gas ap- 
pears to show nearly scale invariant behavior over a very 
broad range of parameters without the need for fine- 



tuning. Understanding this may give insight into re- 
lated problems such as why, in some quantum field the- 
ories with conformal invariance broken by a mass term, 
the sound speed and bulk viscosity remain close to their 
conformal- limit values for a wide range of energies [11]. 

We emphasize that this question is distinct from that 
of small deviations from scale invariance in weakly inter- 
acting 2c? Bosc gases. Quantum gases with an unregu- 
larized delta- function interaction have an SO (2,1) sym- 
metry [12] and exhibit scale invariance. However, the 
cutoff essential to describe an actual short-range interac- 
tion leads to a violation of scale invariance (analogous to 
an anomaly in quantum field theory) and an interaction- 
dependent shift in the breathing mode frequency from 
2ujq [13] in a 2d Bosc gas. The 2d Bose gas experiments 
that sec nearly scale- invariant behavior arc in the weakly- 
interacting regime [14-16], where deviations are expected 
to be small. In contrast, the 2c? Fermi case that we fo- 
cus on is not weakly interacting and we must take into 
account strong interactions. 

Results — We begin by summarizing our approach and 
main results. We consider a dilute Fermi gas in d = 2, 3 
dimensions with a short-range s-wave interaction, aris- 
ing from a broad Feshbach resonance between two spin 
species, each with density n/2. The dimcnsionless in- 
teraction g<i is expressed as (73 = — l/kpa^ in 3c? and 
#2 = log(fcpa2) in 2c?; is the s-wave scattering length 
that sets the dimer binding energy £5 = —l/ma 2 d . The 
fermions have mass to, density n ~ k F , and we set K = 1. 

The quantity of central interest in our analysis is 



ld = (1 + 2/d)P - P (dP/dp) s 



(1) 



which is the deviation of the adiabatic compressibility 
p(dP/dp) s from its value (1 + 2/d)P in a scale-invariant 
system, where the pressure P oc p( 1+2 / d ), Here s = S/N 
is the entropy per particle and p = ran the mass density. 

First, we show that jd governs the difference between 
the frequency uj m of the hydrodynamic monopole breath- 
ing mode and 2ujq for a Fermi gas in an isotropic har- 
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FIG. 1: Schematic plot of the bulk viscosity spectral func- 
tion (solid line) and the scaled "vacuum" contribution 
C(,o(uj)/Co (dashed line) given by (6). The smallness of the 
subtracted sum rule (3) for [£(w) — C£o(w)/Co] implies very 
little spectral weight below the dimer binding energy |ef,| for 
excitations that break scale invariance. 



monic trap V eM (r) = mui^r 2 /2. We find 
cj 2 J4uj 2 = 1- y / d d r 7d (r) I J d d rn(r)V cxt (r). (2) 

Second, we show that jd is related to the exact sum rule 
for the bulk viscosity spectral function C(w): 



(2/tt) / doj [((to) - CCoM/Co] 
Jo 



Id- 



(3) 



Here, C is the "contact" [17, 18] with Co(w) and Co the 
bulk viscosity and contact in the zero-density n — > limit 
at T = 0. The subtraction on the left-hand-side removes 
the large-w tail of (see Fig. 1) and the sum rule thus 
measures the availability of low-energy (< |£f,|) spectral 
weight for excitations that break scale invariance. 



The experimental observations of Ref. [10] in 2d Fermi 
gases imply that 72 <C neF- Our goal is to under- 
stand why 72 is so small for a wide range of interaction 
strengths, even though other measures of the departure 
from scale invariance, such as P — 2e/d (where e is the en- 
ergy density), are not small. 7^ strictly vanishes only at 
the unitary point g 3 = in 3d, and in the weak-coupling 
BCS limit g 2 — > 00 in 2d. However, we will argue that 
there is considerable evidence for an anomalously small 
72 across the entire BCS-BEC crossover. Remarkably, 
within mean field theory [19] 72 IF = for all values of 
172 (and only in d = 2). In addition, the available T = 
quantum Monte Carlo (QMC) data in 2d [20] leads to an 
estimate of 72 that is consistent with zero over the en- 
tire crossover, except possibly near g 2 — 0. We reach the 
same conclusion at finite T using a scaling argument, and 
argue that this is due to the logarithmic dependence of 
52 on density. Using the sum rule (3), we will argue that 
a small 72 also gives insight into the negligible viscous 
damping of the monopole mode. 



Monopole breathing mode — The normal mode solu- 
tions of the hydrodynamic equations with frequency ui 
are obtained from the Lagrangian [9] 



£[«] = uj 2 Jdr Po u 2 (r) - Jdr [p Q l {dP/dp\ (Sp) 
+ 2p (dT/dp) s SpSs + po (dT/ds\ (5s) 



(4) 



describing quadratic fluctuations in entropy 5s and den- 
sity dp about their equilibrium values, sq and Pq. The 
displacement field u(r,i) is related to the velocity v by 
du/dt = v. Conservation of density and entropy gives 
Sp = — V • (pqu) and 5s = —u ■ Vso. Eq. (4) is valid in 
both the normal as well as the superfluid phase, where it 
describes first sound with v n = v s [9]. 

We obtain the result (2) for the breathing mode fre- 
quency using the scaling ansatz u(r, t) = « r cos(wi) 
in (4), together with the Maxwell relation (dP/ds) p = 
p 2 ] (dT/dp) s and the equilibrium identities VPo = 
(dP/dp) s Vp + {dP/ds) p Vs = -noVKxt for V oxt = 
mwgr 2 /2, and VT = {dT/dp) s Vp + {dT/ds) p Vs = 
0. The above scaling ansatz provides a rigorous upper 
bound on the mode frequency [9], Generalizing the vari- 
ational ansatz to« = r Yl n =o u n r2n cos(wt), it is easy to 
show that the corrections to (2) are governed by higher 
powers of 7^. Thus, 7^ rigorously determines the devia- 
tion of the monopole frequency uj m from 2wo- 

We next relate jd to the contact C, given by C = 
2irma2(de / 0,2) s m 2d [21, 22] and C = Anma\(d£ / da^) s 
in 3d [18]. We find 72 = -[C+ f(dC/da 2 ) s }/ / tTrm and 
7 3 = — [C + a 3 (dC/daz) s \ /367rma 3 . This makes it clear 
that uj m = 2wo is strictly valid only for a 2 — > 00, the 
BCS limit in 2d, where C — > 0, and at unitarity in 3d, 
where | a.3 1 — > 00. On the other hand, the breathing mode 
frequency (2) is very sensitive to 7^ 7^ in both 2d and 
3d. Using J d d rnV cxt ~ O(Nep), we estimate that a 
value of 7d as small as O.lnep would give rise to a 5% 
shift in u> m . The fact that no such shift is observed [10] 
in 2d indicates that we must understand why 72 <C ucf 
for a wide range of g 2 and T. 

Viscosity sum rules — The bulk viscosity £ is the only 
transport coefficient that damps the scaling flow u cx 
r [23]. To gain insight into why it is small in 2d, we 
derive a new bulk viscosity sum rule that relates 72 to 
the low-energy spectral weight for excitations that break 
scale-invariance symmetry. 

The bulk viscosity spectral function is re- 

lated by a Kubo formula to the transverse xt (q, w) 
and longitudinal Xi(q> CJ ) current correlators: C{oj) = 
lim 9 ^o rn 2 u! [Im^i — (2 — 2/rf)ImxT] / q 2 . Generalizing 
Ref. [4] to arbitrary d, we obtain the exact sum rule 



2 f°° 

- / duQ{uo) = - (2 - 2j d) X T + X L - pc 2 



(5) 



Here, X T(L) = lim g _> ([j-q, [H, jq]})T(L)/q 2 with the 
current j* = £k CT [( 2k + q) x /2m]cl a c k+(l(7 . The sub- 
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script T{L) denotes the transverse (longitudinal) q —> 

1 /2 

limit [24], and c s = (dP/dp) s is the adiabatic sound 
speed. Evaluating the commutators in (5) for an isotropic 
pair potential with range r , we find the 2d result 
(2/tt) / duj((uj) =2e- pc 2 + aC/m + /3ClnA/m. Here, 
A = l/r is an ultraviolet (UV) cutoff, and a, (3 are con- 
stants. In 3d [4], the terms proportional to C are of the 
form aC/ma?, + /3CA/m. 

The key insight that allows us to obtain physical results 
independent of A is that an UV divergence of precisely 
the same form must arise in the two-body problem. The 
sum rule for Co has the same form as above, but with en- 
ergy density and contact replaced by their zero-density, 
T = values, £o and Co, while c s = for n — > 0. The ex- 
act solution Co of the two-body problem can then be used 
to regularize the divergence in the many-body problem. 
The same idea underlies the standard replacement of the 
bare interaction with the two-particle s-wave scattering 
length in the study of dilute gases [25]. 

The T = 0, zero-density limit Co(^) of the viscosity 
spectral function has an exact representation in terms of 
the sum of all particle-particle ladder diagrams with two 
external current vertices. These are the well-known [26] 
Aslamazov-Larkin, Maki-Thompson, and self-energy di- 
agrams. We thus obtain the 2d result [27] 



Cp Q(u-\e b \) 
4mw \ n 2 (u/\e b \ - 1) + tt 2 



(6) 



for lj > 0, and Co(~ u>) = Co(<^)- In 2d, there is a bound 
state for all values of the scattering length [19]. Thus, 
in the zero-density (single dimer) and temperature limit, 
co = e& = — \jma\ and Co = 47r/a2. The absence of 
spectral weight in Co below |et>| is due to the fact that the 
only excitations in this limit involve pair disassociation 
with a gap |e&| at T = 0. 

The UV divergences can now be removed by look- 
ing at the difference between the sum rule for the in- 
teracting many-body system and that for the T = 0, 
n — > limit, scaled by (C/Cq). In 2c?, we find 
(2/tt) / du: [CM-CCoM/Co] =2e - pc 2 - 2Ce /C . Us- 
ing C = 4:irm(P — e) and Eq/Cq = — l/(47rm), we obtain 
(3) in 2d. The same methodology can be also be used to 
obtain corresponding results for the bulk viscosity in 3d 
as well as the shear viscosity in any d [27]. 

Our main focus will be on (3), which quantifies the 
low-energy spectral weight in C(w) with the high-energy 
tail C(o(w)/Co subtracted out; see Fig. 1. However, we 
can also obtain the total spectral weight in C(w): 



2 f°° 

S M =- I duj((u) = 3P-e-pc 2 s = - 



8nm \dg 2 ) s 



(7) 



S 2 d > for all 32 [21], as required by CM > [4]. In 
Fig. 2, we plot Szd as a function of gi = log^FGte) using 
T = QMC data [20] to evaluate the right-hand-side of 
(7). Both 5*2d and its 3d counterpart [4] (inset of 




FIG. 2: Bulk viscosity sum rule: The 2d sum rule S2d a t 
T = in units of ntF plotted as a function of the coupling 
g2 — ln(kpa2). Inset: The corresponding 3d result [4] as a 
function of gz = — l/(kFaz). In both 2d and 3d, gd +oo 
the BCS limit while n d — > — co is the BEC limit. 



Fig. 2 using the QMC data of Rcf. [28]) are < ne F in 
the BCS region > 1, but become significantly larger 
on the BEC side. 

Apparent scale invariance — We now have all the re- 
sults in hand to discuss deviations from scale invariance. 
So far, we have shown that 7^ controls the deviation of 
the monopole w m from 2wo and also governs the avail- 
ability of low-energy bulk viscosity spectral weight ((lu). 
We can intuitively understand the exact relation (3) be- 
tween a C-sum rule and the shift in the mode frequency 
as a Kramers-Kronig transform of the d.c. bulk viscosity 
C(0) that damps the monopole mode. 

What, if anything, is special about 2d that leads to the 
strong experimental signatures [10] of scale invariance? 
We begin by addressing this question at T = and then 
generalize to finite temperatures. The first clue comes 
from mean field theory (MFT), which in 2d has a trans- 
parent solution [19] across the entire T = BCS-BEC 
crossover: e = ucf/2 — n\eb\/2. This leads to P = nep/2 
and thus r y^ lFT = for all couplings 32- Contrast this 
with the 2d MFT result P — s = n\e b \/2, which is very 
small in the BCS regime but very large on the BEC side. 
This is our first hint of something we will see again: 7^ is 
small in part because it does not involve physics on the 
scale of the dimer binding energy, whereas P — e does. 

To understand how quantum fluctuations beyond MFT 
affect the result for 72, we use T = QMC [20]. We find 
that the QMC-derived 72 is vanishingly small in both 
BCS and BEC regimes, and even for 32 ~ 0, 72 ~ 
(within large error bars) as shown in Fig. 3. We also 
see from this figure that the 2d result is quite different 
from the 3d case. The QMC estimate for 73 (using data 
from Ref. [28]), though quite small on the BCS side of 
the crossover, is large in the BEC region in 3d. 

We now show that the difference between 2d and 3d 
is tied to the form of the dimensionless couplings 32 = 
log(fci?a2) and 33 = — l/kwds- In the BCS limit (ga S> 1), 
the equation of state has the form e = (neF/2)[l+A/gd+ 
B/g^ + . . .] in both 2d and 3d. The pcrturbative Hartree 
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FIG. 3: 72 (solid blue line) and 73 (dashed red line) shown in 
units of ne_F, where n and tF are the two and three dimen- 
sional density and Fermi energy, respectively. The error bars 
on 72 are associated with numerical derivatives of QMC data 
(with errors) [20]. The coupling gd is ln(fci?a2) for d = 2 and 
— l/(kFaz) for d = 3. 



plus "Fermi liquid" corrections arc larger than the pairing 
contribution not shown. (The only qualitative difference 
is that A is negative in 3c? but positive in 2d [29].) Both 
7d, calculated using (I), and P—e = (ad/d) (de/dad), are 
small in the BCS limit. In the BEC limit (gd < and 
|fld| 3> 1), we get £ = — n\sb\/2 + . . ., which is the energy 
density of n/2 dimers with perturbative corrections in 
powers of 1/|<7<z|- The key difference between 2d and 3d 
is in the (^-dependence of the binding energy \sb\, which 
~ exp(|g 2 |) m 2d and ~ 1 / 1 £/3 1 2 in 3d. 

To understand the effects of finite temperature, we 
write the pressure and energy density, related by P — 
n(de/dn) s — e, in the scaling forms P = ne F J r (gd, s) and 
e = nep£(gd, s). There is a qualitative difference be- 
tween the (^-dependence of the scaling functions T and 
E in 2d. The pressure does not have a contribution on 
the scale of the dimer binding energy \sb\ = 1/ma 2 ; i.e., 
it does not have a potentially exponentially large contri- 
bution in ij2 = \og(kpa,2), while the energy density does. 
We have already seen this in the T = MFT results, and 
the same is also observed in the 2d virial expansion [30] . 
We conjecture that the scaling function J 7 is a slowly 
varying function of gi at all temperatures in 2d (except 
in the immediate vicinity of a weak singularity at T c ). 
The equation of state is then P ~ n 2 up to logarithmic 
corrections, leading to a small 72. 

The absence of high energy contributions on the scale 
of the dimer binding energy to P and the compressibility 
is also consistent with 72 being related to the low en- 
ergy spectral weight as shown by our sum rule. Once 
high energy excitations on the scale of \Eb\ are excluded, 
low energy phonons (with a near scale-invariant disper- 
sion Wq(n) ~ y/nq), for instance, dominate the equation 
of state leading to P ~ n 2 and a small 72. Unlike 72, 
however, P — e is not small, as it involves high energy 
contributions on the scale of |ef,| in the BEC regime. 

Another way to characterize the deviation from 
scale invariance, analogous to the "trace anomaly" in 



quantum field theory, is to rewrite (1) as as 7^ = 
-{dP/dg d )P{gd)/d, where f3(g d ) = k F (dg d /dk F ) de- 
scribes the scaling of the coupling gd with respect to 
the momentum scale k F - We see that fi(g-i) = 1 while 
P(93) = 33; reflecting the difference between the loga- 
rithmic and power-law dependence on the density in 2d 
and 3d respectively. In both the BCS (gd 3> 1) and BEC 
(gd <C —1) regions, the 2d beta function is much smaller 
than its 3d counterpart. 

Finally, using the sum rules (7) and (3), we discuss 
the damping of the monopole mode, controlled by £(0). 
Although a small value for the sum rule by itself does 
not rigorously upper-bound C, any physically reasonable 
functional form for the spectral function (e.g., a Drude 
form for w < |e(,|; see Fig. 1) would lead to a very small 
value for £(0). We see from Fig. 2 that in the BCS regime 
9d <] 1, the sum rule J dw£(w) < ntf in both 2d and in 
3d. We would thus expect a very small £ <C n here in 
both 2d and 3d. This sum rule is quite large on the BEC 
side and does not lead to any restriction on £. From 
the low-energy sum rule (3), however, we see that the 
large value of S^d in the BEC limit is entirely dominated 
by the high-energy tail on scales larger than the dimer 
binding energy. Once this is subtracted out, the low- 
energy integrated spectral weight, equal to 72, is very 
small even in the BEC regime (see Fig. I). Thus in 2d, 
we expect the bulk viscosity £ to be very small both in 
the BCS and in the BEC regimes. 

Conclusions — We have shown that the parameter 7^ 
controls the deviation of the breathing mode frequency 
oj m from its scale-invariant value 2wo and also quantifies 
the low-energy spectral weight for excitations that break 
scale invariance, using an exact sum rule. We argue that 
2d is special, with a coupling that depends logarithmi- 
cally on density, leading to a very small 72, even in the 
BEC regime where scale invariance is strongly broken by 
the large dimer binding energy (and hence P—e is large). 
The small 72 also implies, via the 2d sum rule, weak 
damping of the monopole mode in both the BCS and 
BEC regimes. The regime very near gi = deserves fur- 
ther theoretical and experimental investigation, but the 
available evidence suggests that 72 might be small there 
as well. 
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Note added in proof — While completing this 
manuscript we became aware of Ref. [31], which 
analyzes the experiment of Ref. [10] from the different 
perspective of quantum anomalies at T = 0. In the only 
area of substantial overlap, our general result (2) for the 
breathing mode reduces to that of Ref. [31] if we assume 
a polytropic equation of state. 
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FIG. 4: Exact zero-density limit of the current correlation 
function in 2d. Clockwise from top- left: Aslamazov-Larkin, 
the two self-energy, and Maki-Thompson diagrams. Lines 
with arrows denote the vacuum single-particle Green's func- 
tion G; wavy lines, the pair propagator T. The external ver- 
tices are the bare current vertices jx. 



SUPPLEMENTAL MATERIAL 

Viscosity in the zero-density limit — In the zero-density 
limit, the current correlation function can be calculated 
without approximation by summing up all ladder dia- 
grams. (A related approximate calculation for the fi- 
nite density viscosity spectral functions was carried out 
in Ref. [1].) It receives contributions from Aslamazov- 
Larkin (AL), Maki-Thompson (MT), and self-energy 
(SE) diagrams, as shown in Fig. 4 (as long as the en- 
ergy of two particles in vacuum is nonzero and nega- 
tive, as expected in two dimensions [2], the bare loop 
contribution vanishes). Here, the bold lines are zero- 
density Green's functions, G(k) = (ik n — ^k) _1 , with 
£k = £k - M> £ k = k 2 /2m, and fi = -\e b \/2 + 0+. 
5q = Ek<r7z( k > k + q)c k<T Ck+ qC r is the x-component of 
the current operator, with 7 x (k, k') = (k + k%/2m the 
bare current vertex. The internal wavy lines are the vac- 
uum pair propagator T, the inverse of which is given by 



r- 1 (< z ) = -i + i^G(-fc)G(fc 



9) 



(8) 



Using the renormalization scheme [2] 

.g- 1 =T- 1 ( J B)-^(2 ek -2 J B)- 1 , (9) 



where E is some arbitrary (nonzero) energy and 

r- 1 H = £[in(| £b |/c) + H, 



(io) 



in (8) gives the convergent result (after analytic contin- 
uation of the Bose frequency iq m — > uj + i0 + ) 

r-^q, w) = ^ In [1 - (u + i0+ - q 2 /4m + <J/i)/M] ■ 

47T 

(11) 

Here, we have defined fx = —\eb\/2 + 5/j,. Sfj, approaches 
zero from above in the zero-density limit. 

Considering the correlation function for the x- 
componcnt j x of the current, the contributions shown in 



Fig. 4 are given explicitly below. The Aslamazov-Larkin 
contribution is 

Xal(<z) = ~ Yl S *(P> V) T (P) T (P + «)> (I 2 ) 



with the block S x (p, q) of three Green's functions given 
by 

s *m = |E^( k ' k +q) G ( fc +«) G ( fc ) G ^- fc )- ( 13 ) 

P k 

The two self-energy contributions are 

xs(«) = -4 E ^ k + q)r(p)G 2 (fc)G( P -k)x 

p,k 

G(k + q) + (q^ -q) (14) 
and the Maki-Thompson contribution is 

XMT{q) = -■^^7 a; (k,k + q)7 2: (p - k,p- k - q)x 

T( P )G(k)G{ P - k)G{k + q)G{ P -k-q). (15) 

q = (q, iqm), V = (p> ipm), and k = (k, ik n ) are 4- vectors, 
where q m , p m and k n are Bose and Fermi Matsubara fre- 
quencies, respectively. The pre-factors of 2 appearing in 
(13), (14), and (15) arise since there are two spin degrees 
of freedom. 

Using (11), the contributions to the density correlation 
are found after some straightforward, if modestly labori- 
ous, algebra. To order q 2 , the transverse and longitudinal 
components of the MT + SE contributions are 

X T (q^)lMT+SE = 8m ( mcJ )2 9 (^ - N) - (W "> -«), 

(16) 

Xi mt+se = ^—7 xT®^- \e b \) - (u -> -w). 

(17) 

As in 3d, the AL contribution is purely longitudinal [1]: 



xz,(q,w)| AL 



C Q q 2 0(lu - \e b \ 



4m(mw) 2 g(oj) + tt 2 
[l-(g(oj) + 7r 2 )f(iu)]-(LO^-oj). (18) 

In the above, we have defined f(co) = (1 — \eb\/oj) 2 and 
g(uj) = ln 2 (cj/|e(,| — 1). Go = 4:7rm\s b | is the zero-density 
limit of the contact. Using these in the appropriate Kubo 
formulae (see main text, above (5)), one finds (6). Ap- 
plying this procedure to the 3d case, one reproduces the 
C(w) sum rule reported in Ref. [3]. 

Shear Viscosity — In the paper we focus on the bulk 
viscosity sum rule, since this £ is relevant for the damping 
of the monopole breathing mode and has direct bearing 
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on the question of scale invariance. For completeness, we 
briefly discuss here the shear viscosity sum rule. 

We first use the above results to evaluate the zero- 
density, T = limit t] (uj) of the shear viscosity r)(w) = 
lirrig_ > o m 2 wlm.XT / q 2 , with (2 / ir)j dwq{uj) =Xt- For the 
even spectral function [r)o(— w) = r]o(u)], we find 

?7o (w > 0) = (Cb/8mw)(l - \e b \/Lu) 2 e(io - \s b \), (19) 

consistent with the high-frequency spectral tail calcu- 
lated in Rcf. [4]. 

We next write down the convergent subtracted sum 
rule for the shear viscosity. Evaluating the commuta- 
tors in (5) in the main text, one finds (2/tt) J dLurj^uj) = 
e + a'C/m + j3'C In A/m for the shear viscosity sum rule. 
As with the bulk viscosity sum rule, A = 1 /r is an ultra- 
violet (UV) cutoff, and a',/3' are constants. Subtracting 
the sum rule for Ct]q(uj)/Co from the sum rule for the 
full many-body viscosity, one finds 



Following the same procedure to obtain a subtracted 
shear viscosity sum rule in 3d gives the sum rule obtained 
in Ref. [1]. 
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(2tt) / du [r)(u)-Cr) (cj)/C ] =e - Ce /C 
Jo 

=e + C/Airm. (20) 



